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The serotonin 5HT7 receptor has been implicated in numerous physiological and pathological processes from 
circadian rhythms [1] to depression and schizophrenia. Clonal cell lines heterologously expressing recombi- 
nant receptors offer good models for understanding drug-receptor interactions and development of quantitative 
structure-activity relationships (QSAR). Comparative Molecular Field Analysis (CoMFA) is an important modern 
QSAR procedure that relates the steric and electrostatic fields of a set of aligned compounds to affinity. Here, 
we utilized CoMFA to predict affinity for a number of high-affinity ligands at the recombinant guinea pig 5HT7 
receptor. Using R-lisuride as the template, a final CoMFA model was derived using procedures similar to those of 
our recent papers [2, 3, 4] The final cross-validated model accounted for >85% of the variance in the compound 
affinity data, while the final non-cross validated model accounted for >99% of the variance. Model evaluation was 
done using cross-validation methods with groups of 5 ligands. Twenty cross-validation runs yielded an average 
predictive r(q 2 ) of 0.779 ± 0.015 (range: 0.669-0.867). Furthermore, 3D-chemical database search queries de- 
rived from the model yielded hit lists of promising agents with high structural similarity to the template. Together, 
these results suggest a possible basis for high-affinity drug action at 5HT 7 receptors. 

Abbreviations: 5-Carboxarnidotiyptamine (5-CT); 5-HT (5-hydroxytryptamine, serotonin); 5-methoxytryptamine 
(5-MET); 8-OH-DPAT, R- and S-8-hydroxy-dipropyl-amino-tetraline hydro-bromide; ACD, Available Chemical 
Directory; APO, R-apomorphine; bulbocapnine, (-f-)-buibocapnine; CoMFA, comparative molecular field analysis; 
CoMSIA, Comparative Molecular Similarity Indices Analysis; K L , low-affinity agonist dissociation constant, 
CMC, Comprehensive Medicinal Chemistry; lisuride, R-lisuride; mCPP, l-(3-chloro-phenyl)piperazine; May- 
bridge; MDDR, Modem Drug Data Report; NCI, national cancer institute; NPA, R-N-n-propylnorapomorphine; 
QSAR, quantitative structure-activity relationship; SEE, standard error of the estimate; SEP, standard error of 



recombinant receptors 



Summary 



prediction. 



The serotonin-7 (5HT 7 ) receptor represents the most 
recent [1,5] and certainly one of the more intriguing 



Introduction 
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members of the serotonin receptor family [6, 7]. There 
are several reasons for such an interest in the 5HT7 re- 
ceptor. First, it is expressed in many species, resulting 
in cloning from rat brain [1, 8, 9], guinea pig [10], and 
Drosophila [11]. Second, the 5HT 7 receptor appears 
implicated in a number of diseases and physiologi- 
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Tabic J. Affinities of training set compounds at recombinant SHTt 
receptoj-s 



Drug 



Measured K u (M) Predicted Kj (M) 



I mini no c»»t 










S-Mclhn XVtrunt'imino 


_>.U I Li 


— 10 


4.96E 


-10 


5CT 


5.25E- 


— 10 


5.2 I E 


-10 


»'.ri J 


8.7IE- 


- 10 


8.84E- 


-10 


AMnesergicic 


7.80E- 


-08 


Z87E- 


-08 


/a |.»o i no r pn i n e 


5.62E- 


-07 


5.52E- 


-07 


Bromocriptine 


L07E- 


-08 


1 .08E- 


-08 


Bulbocapninc 


2.69E- 


-07 


2.7IE- 


-07 


Ergotamine 


5.0 IE- 


-08 


5.00E- 


-08 


Laurolistine 


3.3 IE- 


-07 


3.34E- 


-07 


R-Lisuride 


9.50E- 


-09 


9.58E- 


-09 


LY2 15840 


1.47E- 


-08 


1.48E- 


-08 


LY53857 


1.02E- 


-07 


9.89E- 


-08 


NPA 


6.92E- 


■07 


6.99E- 


-07 


Oxymetazolinc 


S.01E- 


-08 


4.97E- 


-08 


Sergolcxole 


J.01E- 


-07 


I.02E- 


•07 


R-SOHDPAT 


3.50E- 


•08 


1 .35E- 


•08 


Tcrgtiride 


I.00E- 


■08 


I.63E- 


08 



Table 2. Literature values for affinity of 
8-OH-DPAT 



Isomer 


fCi (nM) 


Reference 


RS 


35 Kj 


9 


RS 


52 Kj 


8 


R 


81-141 Kj 


7 


R 


5623 EC50 cyclase 


7 


R 


5200 EC50 cyclase 


10 


RS 


2300 EC50 cyclase 


1 


RS 


62% Intrinsic activity 


1 


RS 


466 Kj 


11 


R 


39-501 Kj 


6 



the 5HT 7 receptor can proceed more effectively if the 
molecular bases for agonist and antagonist interactions 
with the 5HT 7 receptor are understood. The method 
of Comparative Molecular Field Analysis (CoMFA) 
is one technique that has been shown to be effective 
in increasing the understanding of drug-receptor in- 
teractions at a variety of receptors, [n this procedure 
compounds are aligned according to a specific rule 
(the alignment rule). One approach to this process is 
to align compounds with a template molecule using 
flexible-field lit methods. Then the CoMFA proce- 
dure computes the sleric and electrostatic inlermole- 
cular interaction fields of 3D-grid points surrounding 
the aligned ligands. Finally, the partial least squares 
(PLS) regression method is used to derive an opti- 
mum CoMFA model relating the computed interaction 
potentials to measured, and, subsequently, predicted 
activities. Typically such CoMFA models have high 
predictive validity [2-4]. These models also have use- 
ful qualitative features in that the contour map of the 
steric and electrostatic fields associated with the model 
provides insights regarding the tecepiorsite of interest. 

In the present report we sought to develop a 
CoMFA model for potential agonist compounds inter- 
acting at the 5HT 7 receptor. High-affinity compounds 
having interesting 2D- and 3D- structural features 
were selected and those compounds thought to possess 
at least some agonist activity were retained for a final 
CoMFA. The CoMFA model derived from these com- 
pounds appears to have good cross-validation charac- 
teristics. Furthermore, the model has yielded sufficient 
information to allow formulation of productive search 
queries of several 3D-chemical databases. Finally the 
model shows interesting similarities and contrasts with 
previously derived models for agonist affinity at D, 
and D 2 dopamine 1 2-4] and 5HT, A receptors. 



cal processes including migraine [12], psychosis [13], 
depression [14], gastrointestinal functions [15], and 
eircadian rhythms [I]. The pharmacology of the 5HT 7 
receptor is especially interesting to investigators inter- 
ested in drug development because of its high affinity 
tor a number of antidepressants [14], antipsychotics 
(both typical and atypical) and antiparkinson agents 
[1,8-11]. 

The ability to activate 5HT 7 receptors selectively 
represents a goal with both theoretical and clinical 
significance. Development of compounds to achieve 
this goal and an understanding of the regulation of 



Methods 

Cell culture 

Chinese hamster ovary (CHO) cells iransfected with 
the guinea pig 5HT 7 receptor were grown to near- 
confluence in a modified F-I2 medium, contain- 
ing 10% w/v fetal calf serum, 1% w / v pen j_ 
eill in/streptomycin and 0.3 mg/ml geneticin (5% C0 2 
37 °C) in polystyrene flasks. CHO cells (3-17 x I0 6 ") 
were disrupted in lysis buffer (Tris, MgS04, and 
Na 2 EDTA), membranes were isolated by centriftiea- 
tion at 39 000 g for 1 0 min at 4 °C, and resuspended in 
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CH. — 




-O-CH-CH-CH, 




R-N-n-propylnorapomorphine 
(NPA) 





OCH, 




^J(CH 2 CH 2 CH 3 ) 2 



R - 8OH-0PAT 




I'ifliue /. Ccmlinucd. 



,CH 2 CH 3 



CH*CH 3 



R-TERGURIDE 



STEEDS buffer according to Tsou 1 10]. Membranes 
were frozen in 1 ml aliquots of approximately I mg/ml 
protein concentration. 

Compounds used 

5-Carboxamidotryptamine (5-CT), 5-methoxytrypta- 
mine (5-MeOT), methylsergide maleate, 5-hydroxy- 
tryptamine (5-HT) hydrochloride, and 8-hydroxydi- 
propylaminotetraline hydrobromide(8-OH-DPAT) were 
purchased from Research Biochemtcals International 
(Natick, MA). ( 3 H)-5-CT (50.4 Ci/inmo!) was pur- 
chased from Dupont/New England Nuclear (Boston, 
MA). All chemicaJs used for adenylate cyclase assays 



were purchased from Sigma. Radiolaheled ATP for 
the adenylate cyclase assays was purchased from New 
England Nuclear. 

Compound affinities at recombinant 5HTy receptors 

Affinities of potential ligands were measured at re- 
combinant 5HT 7 receptors as described below |7 1. The 
gene product was originally cloned from cDNA of 
guinea pig hippocampus [10] and stably expressed in 
CHO cells. Some of the compound affinity values have 
been reported previously [7]. Additionally, affinities of 
four additional 5HT 7 receptor ligands were obtained 
from the literature (amesurgide, LY53857, LY2 15840, 
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Displacement of 
3H-CT binding 
to 5HT7 receptors 




Log concentration 



Fi!>inv 2. Displacement of 3H-5CT by 5 ITT ut recombinant 5HT7 receptors in CHO cells. Shown are the results of a single experiment in 
duplicate. 



and sergolexole) [16]. A training set of 17 compounds 
was utilized for the present study. Affinities of these 
compounds are provided in Table 1 . 

For saturation studies using radioligand binding 
methods, membranes were incubated for two hours at 
room temperature with 0.05-10 nM ( 3 H)-5-CT (51.3 
Ci/mmol New England Nuclear, Boston) in STEEDS 
buffer. To each assay tube were added 200 |xl cell 
membrane, 100 jjlI radioligand, 150 |xl buffer and 
50 fx) of buffer (total binding) or 50 |xl of 1 \xM 5- 
HT, to define non-specific binding. For competition 
studies, membranes were incubated with 0. 3-0.7 nM 
( 3 H)-5-CT for two hours with or without competing 
ligands. Duplicate assays were run using ten concen- 
trations of competing ligand at approximately equal 
log intervals from I0" 12 to UP 7 or 10" 10 to \Q- 5 
M depending on the ligand in addition to total and 
nonspecific binding. To each assay tube, 200 |xl cell 
membrane, 100 pJ radioligand, 150 ^tl buffer and an 
additional 50 jxi of buffer (for total binding), 50 [jlI 
of 1 |xM 5-HT (for nonspecific binding) or 50 |xl of 
displacing compounds were added. Assays were ter- 
minated by rapid filtration through 0.3% polyethylen- 
iminepretreated GF/B filters, washed twice with I ml 
ice cold 50 mM Tris-basc, pH 7.4. Bound radioactivity 
was determined by liquid scintillation spectrophotom- 



etry using a Packard Topcount scintillation counter. 
IC50 values were obtained by non-linear regression. 
The inhibition dissociation constant (K,) of each com- 
pound was then determined according to the method 
of Cheng and Prusoff [46]. 

Adenylyl cyclase assay 

Adenylyl cyclase activity was determined according 
to the method of Salomon. [17J There were only 
minor alterations in the protocol. Briefly, CHO cell 
membranes (20 |xg) were incubated for 20 min at 
30 °C in a final volume of 50 |xl of buffer. The 
buffer contained Tris acetate (25 mM, pH 7.6), mag- 
nesium acetate (5 mM), ATP (0.1 mM), cyclic AMP 
(0.05 mM), dithiothrekol (I mM), BSA (0.1 mg/ml), 
GTP (0.01 mM), creatine phosphate (5 mM), creatine 
phosphokinase (50 U/ml), and (ct- 32 P) ATP ( I |xCi). 
The reaction was stopped by the addition of 100 \x\ 
of stopping solution containing 2% sodium lauryl sul- 
fate, ATP (45 mM), and 3'5'-cyclic AMP (1.3 mM). 
( 3 H) Cyclic AMP (~ 10000 cpm) was added to mon- 
itor the recovery of ( 32 P) cyclic AMP. Formed ( 32 P) 
cyclic AMP and ( 3 H) cAMP were isolated by dual- 
column (dowex, alumina) chromatography. Samples 
were boiled for 5 min, each poured over a Dowex 
column that was washed with 0.5 ml HoO, and the 
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A) 



(+} DPAT vs {-) DPAT-INDUCED 
ADENYLYL CYCLASE 
ACTIVITY 



75 



50- 



-25 




mm 
10 -9 



i mnq 

-7 



TTrnaj — 
-6 



"m i rmq — 
-4 



-3 



[DPAT] (M) 



(+) DPAT (n=6): 
EC 50 = O.134nM 
(0.022 - 0.794) 

(-) DPAT (n=5): 
EC 50 = 1.75jiM 
(1.58-2.34) 



B) 



>H v rv. 



R(+)-8-OH-DPAT , 



S(-)-8-OH-DPAT 



Ft nun 
comparison 



(B) stl , lctURll 



effluent was discarded. The Dowex column was then 
washed with 4 ml H>0 into an alumina column, which 
was subsequently washed with 1 .25 ml of O.I M im- 
idazole buffer (pH 7.3), discarding the effluent. The 
column was washed a final time with 2.5 ml imidazole 
buffer, with the effluent collected into a scintillation 
vial with fluor added for counting. The cAMP produc- 
tion was calculated as pmol/mg protein/min. Results 
were expressed as the change in activity, which was 



calculated as ligand-stimulated minus basal activity 
Responses were quantified by dual channel liquid 
scintillation counting using a Packard model 1 900 TR 
liquid scintillation analyzer. Results were converted 
to the mean ± SEM. EC 50 values were calculated 
according to Fleming [18] and reported as the geomet- 
ric mean ±95% confidence intervals. The Student's 
Mest was used to determine statistical differences be- 
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Table Distance information lor training set ligands al recombinant 5HT7 receptors 



Drim N+ - H-bond 1 (N) 


N+ - H-bond 2 


N+ — ccntroid 


N+ 


- plane 


Energy w/ 


energy 












Cheight) 


electrostatics 


(local minimum) 


5CT 5.796 








4.992 


4.35 c-Ot 


26.98 




23.91 


5HT 5.905 








4.884 


9.54 


e-01 


26.64 




23.83 


5MET 5.995 








4.953 


9.96 c-0 1 


31.97 




27.30 


Amesurgidc 5.905 








4.842 


1.08 e+00 


45.57 




38.49 


Aponiorphinc 6.547 mOH 


7.891 pOH 


5.192 


6.61 


e-01 


28.03 




21.35 


Bromocriptine 5.980 








4.888 


3.74 e-01 


39.54 




35.00 


Bulbocapninc 6.482 mOH 


7.833 pOH 


5.173 


1.23 e+00 


39.10 




35.04 


Ergotamine 5.978 








4.891 


4.48 c-01 


60.91 




55.32 


Laurolistine 7.415 mOH 


7.766 pOH 


5.167 


1.39 e-01 


34.60 




16.00 


Lisuride R+ 5.965 








4.883 


5.13 e -01 


33.15 




33.15 


LY2 15840 5.965 








4.883 


5.13 c-01 


57.61 




47.92 


LYS3857 5.963 








4.886 


. 9.03 e-01 


48.10 




41.65 


NPA 6.603 mOH 


7.933 pOH 


5.223 


4.72 c-01 


37.63 




26.44 


Oxymeiuzolinc 7.197 








4.929 


4.64 e-01 


21.10 




19.87 


Sergolexolc 5.979 








4.891 


6.33 e-01 


49.16 




42.08 


R-80HDPAT 4.901 








5.055 


l.38e+00 


38.88 




29.95 


R-Terguride 5.965 








4.886 


7.26e-0l 


33.68 




30.27 


Distances and plane heights are in k. Energy is in kcal/mol. 














Table 4. CoMFA results for iigand affinity at recombinant SHTj receptors 











Fields 


Dielectric 


Min. a 


Grid 


Probe 


#col 


q2 (SEP) 




if 


CoMFA grid 


(E cutofT 


function 




step 


atom type 










box 


kcal per mole) 






size (A) 














Steric 30/Elec 30 


1/r 


2.0 


2.0 


Gsp3+ 


358 


0.583 (1.863) 


6 


Auto 11 


Slcric 30 


1/r 


2.0 


2.0 


Csp3+ 


226 


0.541(1.995) 


4 


Auto 


Elec. 30 


I/r 


2.0 


2.0 


Csp3+ 


202 


0.734(1.488) 


6 


Auto 


Both 20/20 


' l/r 


2.0 


2.0 


Csp3+ 


330 


0.656 (1.691) 


6 


Auto 


Both 10/10 


I/r 


2.0 


2.0 


Csp3+ 


247 


0.734 (1.487) 


6 


Auto 


Both 25/25 


I/r 


2.0 


2.0 


Csp3+ 


351 


0.588 (1.852) 


6 


Auto 


Both 15/15 


I/r 


2.0 


2.0 


Csp3+ 


183 


0.743 (1.463) 


6 


Auto 


Steric 20/Elec 15 


I/r 


2.0 


2.0 


Csp3+ 


307 


0.698 (1.584) 


6 


Auto 


Steric !5/Elec 10 


l/r 


2.0 


2.0 


Csp3+ 


253 


0.745 (1.389) 


5 


Auio 


Steric 15/Elcc 10 


1/r 


1.75 


2.0 


Csp3+ 


274 


0.851 (1.061) 


c 


New region 11 


Slcric I5/Elec 10 


I/r 


1-1.75 


2.2 


Csp3+ 


221 


0.753 (I, 


,434) 


6 


New region 


Slcric 15/Elec 10 


I/r 


.1.75 


1.8 


Csp3+ 


400 


0.833 (1 


.180) 


6 


New region 


Steric 15 


I/r 


1.75 


2.0 


Csp3+ 


223 


0.782 (1.346) 


6 


New region 


Elcc 10 


1/r 


.1 .75 


2.0 


Csp3+ 


165 


0.756 (1.424) 


6 


New region 


Steric 15/Elcc 10 


1/r 


1.75 


2.0 


0 2 - 


277 


0.842 ( 1 


.146) 


6 


New region 


Steric 15/Elcc 10 


1/r 


1.75 


2.0 


H+ 


282 


0,842 (1.145) 


6 


New region 


Steric 15/Elcc 10 


C 


1.75 


2.0 


Csp3+ 


223 


0.782 ( I 


.346) 


6 


New region 


Region Focus 0 


1/r 


' 1.75 


2.0 


Csp3+ 


112 


0.808 (1.265) 


6 


New region 


Sleri 15/Elec 10 




















CoMSlA d 


I/r 


1.75 


2.0 


Csp3+ 


1021 


0.856(1.096) 


6 


New region 



"The auto CoMFA grid box had coordinates X (- 13. 19, 8.84). Y (-10.28, 16.30), and Z ( - 12.79. 5.45). 
h Optimal grid box position was as follows: X (-13.79, 8.54). Y (10.73. 16. 15), and Z (-12.44, 5.45) 
c Defnull discriminant function setting (power = 0.3). 

d CoMSIA included electrostatic, hydrophobic, steric. donor and acceptor components. 
PC = principal components. 
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B) 




tween log EC 50 values for the DPAT (+) and (-) 
enantiomers. 



Computational chemistry 

Overview of the approach 

For all CoMFA, Region Focusing, and CoMSIA 
(Comparative Molecular Similarity Indices Analysis) 
studies, initial conformations for training set lig- 
and alignment were established using CONCORD 
[19, 20]. Conformations of each compound were 



minimized using the MAXIM1N2 procedure within 
[21]. Flexible-field fit methods were then uti- 
lized to determine the conformations of the training set 
compounds that matched the template [2-4]. CoMFA 

c5 S nw, rf0rmed US '" S the Procedures provided within 
oYbYL [ 2-4 1 . 

For the CoMFA we utilized R-lisuride, a high- 
al hnity and conformationally constrained 5HT 7 recep- 
tor agonist, as the template (Figure 1 ). 5-CT also has 
very high affinity for the 5HT 7 receptor. However its 
structure (Figure I) is highly flexible, with a large 
number of possible candidate binding conformations 
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«r W ,v 5 Ali-nniem of compounds in the training set. (A) Stereo view of .U compounds in the training scl. (B) Stereo view of catechol -like 
compounds. (Cj Stereo view of crgoline compounds. (D) Stereo view of indole-hkc compounds. 



possible with a reasonable energy. Therefore, 5-CT 
was less suitable as a CoMFA template than lisuride. 

Distance information 

A crucial aspect of drug binding to serotonin recep- 
tors is the formation of the electrostatic interaction 
between an Asp residue (probably in transmembrane 
region HI of the receptor) and the positively charged 



nitrogen of the iigand [22]. A second essential feature 
of agonist binding consists of the interactions between 
hydroxy I- or other hydrogen bonding elements of the 
compound and various Ser or other residues (likely 
within transmembrane region V of the receptor) [23]. 
Serine residues appear to be particularly useful in the 
monoamine receptors as hydrogen donors or acceptors 
[22, 23]. Important features of an alignment rule in- 
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D) 




Figures. Continued. 



elude Ihe distances from the cationic nitrogen to the 
nitrogen {or other hydrogen binding moiety) of the in- 
dole ring, etc. An alignment rule was obtained from 
those structures having high affinity at the receptor 
(see Table 3; Figure 4). As previously discussed [2- 
4J such general features imply that certain portions 
of an agonist must be capable of electrostatic inter- 
actions and able to form hydrogen bonds to yield both 
high affinity and relative efficacy. CoMFA fields rep- 
resent a portion of the enthalpic aspect of the overall 



free energy of the drug-receptor interaction. Also, the 
relat.onship between free energy and the cquilibrium- 
bmdmg constant is logarithmic. Thus, the affinity 
values given in |xM were expressed as In (l/Ki). 

Initial conformations and template selection 

In the present study, SYBYL [21] v 6.6 was used 
lor the modeling. Initial structures were generated as 
d.scussed above. As previously discussed [2-4] com- 
pounds are unlikely to bind to any receptor in their 
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13 Jj U 15 16 17 18 19 20 tT~ 22 23 

Measured affinity as Ln( 1 /K ) 

Figure 6. Correlation bcLwcen measured and predicted affinities for training set compounds. 



global minimum energy conformations. The reason 
is that some extent of torsional change or flexion of 
rotaiable bonds is required to adapt the ligand and re- 
ceptor to a complex of lower energy. Therefore, the 
so-called 'minimum' energy conformation resulting 
from a MAXIMIN2 procedure is only a useful starting 
point for possible candidate (binding) conformations 
of the compound. However, it is important to restrict 
the possible ligand conformations to those that are en- 
ergetically accessible. Generally, a 10 kcal/mol cutoff 
value (difference belween energies of the minimized 
and aligned conformers) is considered reasonable in 
CoMFA studies [24, 25]. Thus, energies of bound 
conformations did not exceed 10 kcal/mol above the 
'minimum' calculated by the MAX1MIN2 procedure 
within SYBYL. 

Alignment 

Flexible field-fit protocols represent an improvement 
over 'static* lit methods for compound alignment. This 
is because they allow minor changes in compound 
conformation as a result of the field fit process 126]. 
Flexible field fits were done using the MAXIMIN2 
procedure within Sybyl with the region set to that 
for the default CoMFA (for initial fits) and to the re- 
gion associated with the best CoMFA (for subsequent 
fits). Steric and electrostatic fields extracted for the 
template compound (lisuride) were used in the fitting. 
AMI charges were recalculated following the flexi- 
ble field fit using the MOPAC procedure and the new 



conformational energy of the compound was obtained. 
The new charges were assigned and the compound 
replaced in the molecular data base. Following each 
flexible field-fit, the CoMFA model was derived using 
the SYBYL default parameters (see below). 

Alternate compounds alignments 

Amesergide was realigned using the flexible field fit 
procedures described above (Figures 7A and 7B). The 
energy of the compound was 45.3 kcal/mol (E m i n = 
42.3 kcal/mol). The alternate version of amesergide 
was substituted for the original version in the molecu- 
lar spread sheet and autofill performed for the CoMFA. 
The q 2 was reduced slightly (from 0.851, shown in 
Table 4, to 0.833). 

Similarly, laurolitsine was reoriented to make the 
other pair of hydroxy 1 groups available for hydrogen 
bonding to the receptor (Figures 7C and 7D), aligned 
with the template compound and substituted for the 
original version of laurolitsine. The energy of the 
compound was 15.8 kcal/mol (Emin. 7.6 kcal/mol). 
Substitution of the new version of laurolitsine and re- 
calculation of the CoMFA/PLS lowered the q 2 to 0.705 
(relative to the original 0.851 ). 

CoMFA and partial least squares (PLS) 

We have systematically investigated the effects of 
changing several CoMFA parameters, including di- 
electric (function of 4 l/r' vs. constant), grid step 
size (1.5, 2.0, 2.5 A), probe atom type (H + , 0~, 
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C^), and the column filtering values (from 1.0 to 
4.0 kciil/mole in 0.5 kcal/mole increments). The initial 
cross-validation was performed using the 'leave one 
out* method. 

Residual values for In (i/Kj) represent the differ- 
ence between the affinity predicted by the model vs. 
that measured. Training set compounds having high 
residuals were field-fit again and the CoMFA again 
determined. As described this is a 'manual-iterative' 
procedure in which the field- (it is refined based on 

'in CoMFA, the probe ulom interne Is with each point in the 3D 
grid box. Al some of those points no atoms can be found. At oth- 
ers, atoms arc associated with training set compounds having low 
affinity. Al still other grid points the probe atom would interact with 
atoms of compounds having high affinity. A variety of probe atom 
types are possible. Tiie telrahcdrnl carbon with a positive charge is 
one of the most commonly used. 



improvement of the predictive validity of the CoMFA 
model until no further improvements are attainable 
[2-4]. 

PLS is a regression method that is useful in sit- 
uations in which the number of predictor variables 
greatly exceeds the number of 'subjects 1 (compounds) 
[27]. As applied to computational chemistry, PLS pro- 
vides an extremely useful strategy for inclusion of 
detailed information about the charge and bulk char- 
acteristics of biological molecules [26, 28-32]. PLS 
routinely incorporates thousands of predictor variables 
(as the steric and electrostatic interactions of a probe 
atom with each point in a 3D grid box) for a train- 
ing set of compounds. Predictors are determined at the 
grid points of a 3-dimensionaI lattice (box) containing 
all members of the training set of aligned compounds. 
As in all regression methods, the PLS retains only 
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those predictors in the final model that substantially 
account for the variance of the dependent variable 
(drug affinity). For example, the cationic nitrogen is 
clearly a necessary feature of drug binding to sero- 
tonin receptors. However, this feature has essentially 
no weight in the final CoMFA model because it is 
common to all training set compounds. 



Region focusing, and CoMSIA 

We tested the CoMFA using the Region Focus fea- 
ture within SYBYL with the default SYBYL settings 
(discriminant function to the 0.3 power). This method 
uses a discriminant (or alternate) function to select 
those CoMFA components from an optimal non-cross- 
validated PLS analysis that contribute substantially to 
affinity prediction. Computational chemistry proce- 
dures for determining drug-receptor interactions only 
analyze a portion of the enthalpic aspects of bind- 
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ing. Other contributors to ligand binding, such as 
solvation, are ignored. However, a convenient first ap- 
proximation of a portion of the solvation interaction 
may be done by incorporating hydrophobicity into the 
CoMFA model. We examined the effect of the hy- 
drophobicity variable on the final model, by means of 
a variant of the standard hydrophobicity term derived 
from CoMSIA [4]. This feature was incorporated into 
the molecular spreadsheet. 

Alternate template 

We also tested the rigor of the CoMFA model by 
using a highly flexible compound as the template. 5- 
MET has the highest affinity of any compound in the 
training set (Table I). However, it is highly flexible 
(Figure 8). For this reason, it would not normally be 
considered a likely candidate for a CoMFA template 
1 32]. However, we chose the local minimum con- 
centration for this compound. All other training set 
compounds were aligned with this alternate template 
and the CoMFA optimized as described above. 

Prediction of affinities for a subset of training set 
compounds from a CoMFA model derived from the 
remaining compounds 

Affinities for four of the training set compounds (ame- 
sergide, sergolexolc, LY2I5840, and LY 53857) were 
derived from the literature [ 1 6]. Therefore, it was use- 
ful to attempt to demonstrate comparability of affinity 



measurements across laboratories by predicting affini- 
ties of these four compounds from a CoMFA model 
based on the remaining 13 compounds of the training 
set. Affinities were predicted using the predict prop- 
erty feature of the Sybyl molecular spreadsheet from 
the non-cross validated PLS analysis in two ways. 
First, we performed the PLS using only the 13 desig- 
nated compounds but the original optimized CoMFA 
model. Second, we derived a new CoMFA/PLS analy- 
sis based entirely on the 13 remaining compounds. 

3D-Search of chemical databases 

One useful way to evaluate alignment rules associ- 
ated with CoMFAs is to perform flexible 3D searches 
of large chemical databases for potential Mead' com- 
pounds that possess the desired features [33 J. Here, we 
used the UNITY software accessible through SYBYL 
to perform searches of the following databases: NCI, 
MDDR, ACD, CMC, and Maybridge [3, 20]. The 
query was a simple alignment rule representing dis- 
tance between the cationic nitrogen, the hydrogen 
bonding group mela- to that nitrogen, and the C2 car- 
bon of the aromatic ring. Distance constraints were 
based on the highest affinity ergoline compounds in 
the training set for these searches [3]. Subsequently, 
results from the searches (UNITY hit lists) were 
loaded into SYBYL molecular spreadsheets and a 
SYBYL similarity search performed with R-lisuride 
as the search template. This allowed a convenient rank 
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Figure 9 Contour map relating compound nfflniiy ai5HT 7 receptors to electrostatic and steric intcrmolccular interaction fields using n lisuride 
template. Positive electrostatic charge is favored (blue) or not favored (red) for high-affinity (80:20). Steric bulk is favored (green) or not favored 
(yellow) for high-affinity (80:20). 



ordering of the compounds in the hit list based on their 
structural resemblance to lisuride. 

Results 

Drug affinity 

We determined the affinity of 13 of the 17 training 
set compounds at 5HT 7 receptors stably expressed in 
CHO cells as described using duplicate 10 point dis- 
placement curves against (3H)-5CT [7]. Competing 
ligands were run at approximately equal log inter- 
vals over approximately 5 orders of magnitude. Fig- 
ure 2 shows data from a sample experiment for 5HT. 
Assays were performed under conditions that pre- 
vented coupling of the receptors to G proteins. Thus, 
Hill coefficients for competition binding curves did 
not differ significantly from unity (data not shown). 
Affinities of four additional training set compounds 
(amesurgide, LY2 15840, LY53857, sergolexole) were 
obtained from the literature [16]. Affinities of all 17 
compounds are given in Table I . 

Adenylate cyclase activity 

Figure 3A compares the ability of R- and S-8-OH- 
DPAT to elevate adenylate cyclase activity in mem- 
branes prepared from CHO ceils heterologously ex- 
pressing 5HT 7 receptors. The R-isomer was approx- 
imately an order of magnitude more potent than the 



S-isomer of 8-OH-DPAT. The derived CoMFA model 
was successful in determining the affinity of the ac- 
tive isomer. The more active R-isomer (in terms of 
ability to elevate adenylate cyclase) had "an EC50 
of 130 nM vs. ^1750 nM for the less active S- 
isomer (Figure 3A). Figure 3B compares the structures 
of the active R- and less active S-isomers of 8-OH- 
DPAT. The EC50 found in the present experiments 
is lower (higher potency) than those previously re- 
ported in the literature (Table 2). This may reflect, in 
part, the density of receptor expression or the coupling 
efficiency between receptor and signal transduction 
proteins derived from CHO cells. 

Two-dimensional structures and training set 
compound alignment 

We used flexible field fit procedures to determine 
the final alignments and conformations of training 
set compounds. The cationic nitrogen and hydrogen 
bonding regions were important features of the final 
alignment. 

Alignment rule for SHTj receptor affinity 

Figure 4 shows the alignment rule for binding of lig- 
ands to 5HT 7 receptors. Panel A indicates the typical 
features, including the cationic nitrogen, hydrogen- 
bonding regions, and the centroid of the five- 
membered ring. The minimal height of the cationic 
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nitrogen above the plane of the five-membered ring is 
shown more clearly in Figure 4B of the same figure. 

Table 3 shows the information derived from the 
alignment rule for the compounds of the final train- 
ing set. R-lisuride was the template for the final 
CoMFA because of its structural similarity to other 
high-affinity compounds in the set, its relatively 
high-affinity at 5HT 7 receptors, and its relative lack 
of conformational flexibility. The highest affinity 
compounds (5HT, 5-methoxytryptamine, and 5-CT) 
demonstrated a distance between the cationic nitrogen 
and the carbon of the five-member ring 'meta' to it of 
5.8-6.0 A. The height of the cationic nitrogen above 
the plane of the indole ring is another important align- 
ment feature. The value for this element ranged from 
0.4- 1 .0 A for the high-affi nity compounds at the 5HT 7 
receptor. 

Ligand alignment via flexible field-fit 

Figure 5 presents the alignment of the training set 
hgands. Figure 5A presents a stereo view of the en- 
tire training set. Figure SB shows selected training set 
compounds having a catechol-like structure, Figure 5C 
shows compounds having the ergoline structure, while 
Figure 5D shows indole compounds. Taken together 
these four panels demonstrate the significant similarity 
of key binding features obtainable for compounds with 
extensive 2D-structural diversity. 

Comparative Molecular Field Analysis/Partial Least 
Squares regression (PLS) 

PLS incorporates as potential predictor variables the 
steric and electrostatic interactions of a probe atom 
at the grid points of a 3-dimensional lattice (box) 
containing all members of the training set of aligned 
compounds. PLS retains only those predictors in the 
final model that substantially account for the variance 
of the dependent variable (drug affinity). PLS methods 
specifically deal with common situations in compu- 
tational chemistry in which there are thousands of 
predictor variables for a few compounds [28-3 1 ] 

The final CoMFA model used both steric and 
electrostatic features. A steric cutoff of 15 and an 
electrostatic cutoff of 10 kcal/mol were used in the 
final model. The distance-dependent dielectric func- 
tion yielded an improvement over models holding 
this value constant. The final CoMFA model utilized 
a Csp3-|- probe atom, and a 2.0 A step size. The 
SYBYL default region (within the grid box) was al- 
tered as described in Table 4 to improve the CoMFA 



Table 5. Prediction of aflinilies of amcsergidc sergolex- 
ole, LY2I5840, and LY53857 from the CoMFA model 



Compound 


Properly 


Value 


CoMFA model 


Amesurgide 


Predicted 


16.41 


Optimized 11 




Measured 


16.37 


Optimized 


Sergolexolc 


Predicted 


16.56 


Optimized 




Measured 


16.11 


Optimized 


LY2 1 5840 


Predicted 


16.47 


Optimized 




Measured 


18.04 


Oniirniviv I 


LY53857 


Predicted 


16.92 


Optimized 




Measured 


16.09 


Optimized 


Amesurgide 


Predicted 


16.11 


Default 1 * 




Measured 


16.37 


Default 


Scrgolexole 


Predicted 


16.15 


Default 




Measured 


16.11 


Default 


LY215840 


Predicted 


15.61 


Default 




Measurecl 


18.04 


Default 


LY53857 


Predicted 


16.40 


Default 




Measured 


16.09 


Default 



. „,,. iiiuuci irom lame 4 with 

that CoMFA region. The PLS for 1 3 compounds remain- 
ing in Hie training scl yielded a q 2 of 0.904 (5 principal 
components. SEP = 1.035, and 246/3024 columns used 
in the analysis). 

Predicted using a new CoMFA model based on the re- 
maining 13 compounds of the training set. The PLS for 
13 compounds remaining in the training set yielded a 
q- of 0.656 (5 principal components, SEP = 1.959, mid 
325/2772 columns used in the analysis). 

model. The minimum a value in the final model was 
1.75 kcal/mole. Together, these parameters yielded 
an optimal final CoMFA model. Because these para- 
meters were similar to those for previously generated 
CoMFA models for agonist affinity at dopamine re- 
ceptors, they also facilitated comparisons among the 
models (see below) [34]. 

The cross-validated r 2 (q 2 ) values that resulted 
from the various CoMFA options for in (I/Kj) as 
the target property are given in Table 5. The optimal 
CoMFA/PLS run for the 17 compounds of the train- 
ing set using the SYBYL default settings yielded a q 2 
value of 0.85 1 with 5 principal components and a .stan- 
dard error of prediction (SEP) of 1 .061. The PLS with 
5 components using the no-validation option yielded 
an r- of 0.996 and a standard error of the mean (SEM) 
at 0.183 with steric features of the model accounting 
tor 4 1 % and electrostatic features accounting for 59% 
of the variance in the dependent variable. 

Region Focusing within SYBYL was conducted on 
the optimal non-cross validated CoMFA model This 
protocol includes only those features of the CoMFA 
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that optimally contribute to the model. As expected, 
*an examination of Table 4 shows that the number of 
CoMFA columns used in the model resulting from Re- 
gion Focusing was somewhat smaller than the number 
of columns used in the original optimal model (112 
of 3024 vs. 274 of 3024). Region focusing (using the 
default setting) slightly worsened the q 2 to 0.808 (with 
6 principal components and a SEP of 1 .265). 

Figure 6 shows the relationships between mea- 
sured and predicted affinities for the optimum CoMFA 
model. Two of the training set compounds had pre- 
dicted affinities outside the 95% confidence intervals. 
They were 5HT and 5 MET. An initial estimate of 
a possible entropic contribution to drug binding was 
made using a hydrophobic term added as a component 
of a CoMSIA along with the steric and electrostatic 
and donor and acceptor atoms (within SYBYL 6.6). 
This resulted in a PLS in which 5HT 7 affinity was pre- 
dicted by a combination of five terms: electrostatics, 
steric features, hydrophobics, donor atom character- 
istics, and acceptor atom characteristics. As shown 
in Table 4, use of these CoMSIA features (at default 
settings) resulted in a q 2 value of 0.856 (6 principal 
components, an SEP of 1 .096, and use of 1021 out of 
7560 columns in the analysis). 

Alternative compound alignments 

We explored alternative alignments for two com- 
pounds for which the alternative binding modes 
seemed most reasonable. First, we utilized a dif- 
ferent alignment for amesergide to allow the third 
nitrogen to serve as a hydrogen bonding moiety. Fig- 
ure 7A compares the orientations for amesergide in 
the original and alternative alignments. Figure 7B 
shows amesergide in relation to the template com- 
pound (lisuride). Amesergide was aligned using the 
flexible field fit procedures described above. The en- 
ergy of the compound was 45.3 kcal/mol (E m in = 
42.3 kcal/mol). The alternate version of amesergide 
was substituted for the original version in the molecu- 
lar spread sheet and autofill performed for the CoMFA. 
The q 2 was reduced slightly (from 0.851, shown in 
Table 4, to 0.833). 

Similarly, laurolitsine was reoriented to make the 
other pair of hydroxyl groups available for hydrogen 
bonding to the receptor (Figures 7C and 7D), aligned 
with the template compound and substituted for the 
original version of laurolitsine. The energy of the 
compound was 15.8 kcal/mol (E m in, 7.6 kcal/mol). 
Substitution of the new version of laurolitsine and re- 
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calculation of the CoMFA/PLS lowered the q 2 to 0.705 
(relative to the original 0.85 1). 

Alternate template 

We used the high affinity compound 5-MET as an 
alternate template for the CoMFA. We selected the 
local minimum energy conformation as the template 
conformation of this compound. Other CoMFA/PLS 
procedures were done as with the original model using 
lisuride as the template. Using the default Sybyl pa- 
rameters, we obtained a q 2 of 0.804 (with 5 principal 
components, an SEP value of 1 .2 1 6, and incorporation 
of 362 of 3024 columns into the model). Thus, the data 
suggest that it may be useful to use the highest affinity 
compound in the training set as an alternate template 
for compound alignment. 

Prediction of affinities 

Prediction of affinities was done for a subset of train- 
ing set compounds from a CoMFA model derived 
from the remaining compounds. Affinities for four of 
the training set compounds (amesergide, sergolexole, 
LY2 1 5840, and LY 53857) were derived from the liter- 
ature [16]. We evaluated the comparability of affinity 
measurements across laboratories by predicting affini- 
ties of these four compounds from a CoMFA model 
based on the remaining 13 compounds of the train- 
ing set. Affinities were predicted using the predict 
property feature of the SYBYL molecular spread- 
sheet from the non-cross validated PLS analysis in 
two ways. First, we performed the PLS using only 
the 13 designated compounds but the original opti- 
mized CoMFA model. Second, we derived a new 
CoMFA/PLS based entirely on the 1 3 remaining com- 
pounds. These results are summarized in Table 5. 
Using a PLS derived from the original CoMFA model, 
prediction of affinity was extremely good for ame- 
sergide, sergolexole, and LY53857 (<0.5 log unit). 
Prediction was not as good for LY2 15840 (measured 
affinity was 18.04 vs. predicted affinity of 16.47 Ln 
units). Similarly, another set of predictions for affini- 
ties of the same compounds was derived from a new 
CoMFA model based on just the 13 remaining com- 
pounds of the training set. Once again, predictions 
were excellent for 3 of the 4 compounds and less 
good for LY2 15840. Overall, the results indicate that 
measurement of drug affinity across laboratories was 
consistent. 
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Contour map 

Figure 9, shows the contour map of the final 5HT/ 
CoMFA model. Both steric (80:20 as favored: unfa- 
vored) and electrostatic (80:20) fields were utilized, 
since this improved the standard error of prediction 
(SEP) and the q 2 value over alternative models. The 
contour map showed several potentially important 
characteristics. First, it exhibited little in the way of 
electrostatic features in the vicinity of the cationic ni- 
trogen. This is consistent with the idea that differences 
among the training set compounds in this region of the 
molecule are not of major importance in high affinity 
binding. This is consistent with the importance of the 
electrostatic bond between the N+ and an Asp (likely 
within transmembrane region III) in binding of all 
5HT 7 ligands [22]. Steric features near the cationic ni- 
trogen are of some importance in differentiating ligand 
binding affinity, particularly sterically favored regions. 

A second salient characteristic is that hydrogen- 
bonding regions am highly significant predictors of 
high affinity for 5HT 7 receptor ligands. The regions 
in which positive charge is favored include the vicin- 
ity of the hydrogen-bonding nitrogen and the regions 
near both the five- and six-membered rings. Steric fea- 
tures, especially where steric bulk is favored for high 
affinity, appear prominently in the model, particularly 
a sterically favored region that lies close to the cationic 
nitrogen. 

Figure 10A-E, shows the contour map with ac- 
ceptor, donor, electrostatic, steric, and hydrophobicity 
components resulting from a CoMSIA for the training 
set. Here, the PLS was run using the hydropho- 
bic, steric/electrostatic, and donor/acceptor CoMFA 
columns as independent variables with Ln (l/Kd) as 
the dependent variable. All contours are shown as 
favored: unfavored. A simple interpretation of this 
contour map is that the hydrophobically favored re- 
gions represent those portions of the ligand that are 
least exposed to the surface of the receptor during 
binding. A comparison of panels A and B in this 
figure indicates that hydrophobicity complements pos- 
itive electrostatic charge in yielding information about 
likely hydrogen bonding regions of the receptor ac- 
tive site that interact most favorably with the ligand. 
Specifically, it appears that the alkyl extension of the 
lisuride template molecule is associated with an elec- 
trostatically favored region (blue region in the upper 
left portion of Figure IOC). This site may be posi- 
tioned in relation to the active site of the 5HT 7 receptor 
by the presence of the hydrophobically-favored (or- 



Tubh6. Cross-validation analysis results for CoMFA model from 
Tabic 4 al recombinant 5HT7 receptors 



Crass- validation 


1 


# principal 


Sum doitl error of 


run # (5 groups) 




components 


prediction 


I 


0.734 


4 


1.358 


2 . 


0.866 


4 


0.964 


3 


0.669 


4 


1 .5 1 5 


4 


0.741 


3 


1.287 


5 


0.8 IK 


4 


1.122 


6 


0.847 


6 


1.127 


7 


0.707 


4 


1 .426 


8 


0.717 


4 


1 .399 


9 


0.736 


4 


1 .352 


JO 


0,748 


5 


1.380 


1 1 


0.867 


5 


1.002 


12 


0.836 


4 


1.066 


13 


0.738 


|» 


1.348 


14 


0.734 


6 


1.487 


15 


0.724 


4 


I.3H3 


16 


0.814 


5 


1.187 


17 


0.860 


4 


0.984 


18 


0.848 


5 


1.073 


19 


0.854 


4 


1.007 


20 


0.724 


6 


1.5 14 



PLS cross-validations runs were done using the optimum model 
from Table 4. This model had a steric cutoff of 15 and an electro- 
static cutoff of JO kcal/mol. The model also employed t\ minimum 
sigma or 1.75 A, a step size of2.0 A. a dielectric equal to l/r, and 
a Csp3+ prone atom. The model used a CoMFA grid box with 
coordinates X (-13.79, 8.54), Y (-10.58, 16.15). and Z (-12.44, 
5.45) with a total of 3024 possible CoMFA columns. 



ange region in the upper left portion of Figure I0E) 
portion of the molecule. More generally, those re- 
gions in Figure I0E in which hydrophobicity is not 
favored for high-affinity binding (white in Figure 10E) 
serve to define the regions in which positive charge 
favors high-affinity (blue in Figure 10C). Additional 
hydrophobic capacity favored (orange in Figure IDE) 
and not favored (white in Figure I0E) in some of the 
regions in which steric bulk favors high-affinity bind- 
ing (green in Figure 10D) further served to refine the 
steric features of the model. Qualitative (contour map) 
aspects of the CoMFA model for high affinity binding 
at the 5HT 7 receptor appear to be enhanced by the 
addition of a hydrophobic component. 

Evaluation of the CoMFA model 

To evaluate the optimum CoMFA model, a series of 20 
cross-validations was runs using groups of 5 per run. 
Results are summarized in Table 6. The mean q 2 was 
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Fivure JO, CoMSIA contour maps. Acceptor sites: Favored = magenta; unfavored = red (60:40). Donor sites: Favored = cyan: unfavored = 
purple (60:40). Electrostatic: Favored = blue; unfavored = red (60:40). Sieric: Favored = green; unfavored = yellow (60:40). Hydrophobic: 
Favored = orange; unfavored = while (60:40). 
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Figure 10. Continued. 




0.779 ± 0.015 with a range from 0.669 to 0.867. The 
mean SEP for the 20 runs was 1.249 ± 0.043 with a 
range from 0.964 to 1 .5 1 5. Thus, the present model is 
robusi enough lo predict ligand affinity at 5HT 7 recep- 
tors even when using information from only 12 of the 
17 compounds in the training set. 

JD-chcmical data base search and similarity index 

Intially, the validity of the query (Figure 11 A) was 
tested by applying it to the training set. The 'hits' 
from the search included all four of the highest affinity 
compounds of the training set (5-MET, 5-CT, 5HT, 
and lisuride) and excluded four of the lowest affinity 
compounds (apomorphine, bulbocapnine, laurolitsine, 
and NPA). 

Tabic 7 summarizes the results from searching five 
3D-chcmicul databases via the UNITY program acces- 
sible within SYBYL [3] with a simple search query 
derived from the alignment rule for the 5HT 7 recep- 
tor. Figure 7A shows the search query. In all instances 
a significant reduction in the number of compounds 
(Table 6, column 2 vs. column 3) was obtained while 
all database searches yielded at least two hits. Thus, 

2 Thc 3D-searchcs were clone using an earlier version of the CoMFA 
model widi 15 ligands (K-OH-DPAT and tcrguridc were not in- 
cluded). 



Table 7. 3D-chcmical database search results for lig- 
ands at recombinant 5HT7 receptors" 



Database 


Compounds 


Hits 


Similarity 


Maybridge 


61184 


9 


0.49 


ACD 


201540 


51 


0.52 


CMC 


7100 


2 


0.45 


MDDR 


90158 


38 


0.79, 0.74 h 


NCI 


I 17649 


84 


0.50 



a The 3D-searches were based on a preliminary 
CoMFA model without 8-OH-DPAT and tcrguridc, 
" h Seurch of this database yielded compounds #148982 
and// 147760. 



a simple 3D search query yields a reasonable number 
of potential lead compounds for the discovery of novel 
5HT 7 receptor ligands. Refining the search query us- 
ing an exclusion volume derived from the CoMFA 
contour map can be used to further refine the hit lists 
resulting from the search. 

Similarity indices were calculated for the com- 
pounds in each hit list by using the similarity index 
features within the SYBYL molecular spreadsheet 
subroutine (Table 7, column 4). Here, the refer- 
ence compound for the ranking was the template, 
R-lisuride. In most cases, the compounds yielded 
modest similarity to lisuride (in the range of 0.5 on 
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a 0-1.0 scale). However, one compound of inter- 
est was identified as a result of searching the hit 
list derived from the MDDR database. This com- 
pound had a similarity index of 0.79 (compound 
#148982; Figure MB) compared to lisuride. A sec- 
ond compound (#147760) had a similarity index 
of 0.74 relative to lisuride. These two compounds 
are depicted together in Figure |]Q Compound 
#148982 was identified as l-(6-Methyl-9-ergo)en-8- 
beta-y-l-methyI)-imidazolidine-2, 4-dione N6-oxide. 



Compound #147760 was identified as l-(6-Methyl-9- 
ergolen-8-beta-y-l-methyl)-pyrimidiiic-2, 4( 1 H,3H)- 
dione N6-oxide. The source for both compounds was 
Farmitalia Carlo Erba. Both compounds were indi- 
cated to have anti-hypertensive activity. 

Comparison with other CoMFA models 

We have previously reported CoMFA models for G- 
protein coupled dopamine receptors positively (D|, [2, 
3J) and negatively (D 2 , [3, 4]) linked to adenylate 
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cyclase. The models shown were obtained under con- 
ditions virtually identical to that for the present 5HT7 
model. This facilitates a comparison among them. 

Figure 1 2 presents contour maps for Di and D? re- 
ceptors. All are shown with 80:20 (favored unfavored) 
lo facilitate comparisons among them. The current 
5HT 7 contour map (Figure 12A, the same contour as 
Figure 9 but slightly rotated) and the D| map (Fig- 
ure 12B) are similar. Similarities are most marked in 
ihe presence of positive charge being favored near the 
hydrogen bonding regions of the molecule and in the 
general significance of steric aspects of the intermole- 
cular interactions. However, the 5HT 7 contour map 
differs from that for the Di dopamine receptor. For 
the former, steric bulk is highly favored over most 
regions of the receptor for high-affinity drug-receptor 
interactions, whereas for the Di receptor steric bulk is 
generally not favored [2, 3]. 

Comparing the contour maps for the D2 dopamine 
receptor (Figure 1 2C) with that for the 5HT 7 receptor, 
a similarity is that positive charge is highly favored - 
near the hydrogen-bonding region for both models. 
However, the D2 model has a second large region at 
the opposite end of the binding pocket in which posi- 
tive charge is also favored, while the 5HT 7 model does 

3 Thc images used in this figure were generated for Lhc present paper 
from the data used to generate Figures 5 and 6 in Wilcox et al. 1 31 
CoMFA-hascd prediction of agonist affinities at recombinant DI vs. 
D2 dopamine receplors. 



not. Sterically, the two models are dramatically differ- 
ent from one another. Thus, the 5HT 7 model contains 
a large spatial region in which steric bulk is favored 
for high affinity binding while the D2 model generally 
is unfavorable for high-affinity in the same relative 
position. [3, 4] 

Discussion 

Data shown in Figure 3 represent the first report of 
the potency and intrinsic activity of the resolved iso- 
mers of the 5HT receptor agonist, 8-OH-DPAT. As 
indicated in Table 3, the information on the racemate 
suggests that this compound is a potent partial ago- 
nist at 5HT7 receptors as well as 5HTia receptors 
CpKa 8.2) [35]. The fact that the potency is compa- 
rable at 5HT 7 and 5HTia receptors (present results 
and [7]) suggests a role For 5HT 7 receptors in many 
of the behavioral actions of 8-OH-DPAT. One of the 
highest affinity compounds in the training set is 5- 
CT. This drug exhibits high affinity for members of 
the 5HTi subfamily, but markedly lower affinity for 
the 5HT 2 , 5HT}, and 5HT 4 receptor subfamilies [35]. 
These results suggest that the 5HT| A receptor (neg- 
atively coupled to adenylate cyclase) and the 5HT 7 
receptor (positively coupled to the enzyme) share key 
features for agonist binding [34]. Further support for 
the similarity in active sites for the 5HTia and 5HT 7 
receptors is provided below. 
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is lavored (hluc) ornol lavored (red) lur high-ariinity (80:20). Stcric bulk is favored (green) or not favoivd (yellow) ror liieh-nffiniiy (80 ■->()) 
(This ligurc is minted from thai shown in Finnic 9.) (B) D, ilnpuminc receptor model - Br- APR template. R-lisuridc is shown in ilic (mure' 
Positive electrostatic charge is favored (blue) or not favored (red) tor high-affinity (H0:20). Sieric bulk is favored (green) or not favored (yellow) 
lor high-afliniiy (80:2(1). (Cj D 2 dopamine receptor model - Bromocriptine template, R-lisuridc is shown in the ligurc Positive elect rnst-ilic 
charge is favored (bine) or not favored (red) for high-nffinily (80:20). Stcric bulk i.s favored (green) or not favored (yellow) for high-ifliniiy 
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The present results add to the information on in- 
teractions of ligands with G protein coupled receptors 
obtained through computational chemistry methods 
1 2-4, 22, 26, 36 1. Here, we determined the extent to 
which CoMFA-derived models of the 5HT 7 receptor 
could usefully predict affinity for a highly diverse set 
of ligancls. Using a series of 20 cross-validation runs 
with groups of 5 compounds to evaluate the CoMFA 
model we found that the mean q 2 was 0.779 ± 0.015 
with a mean SEP of 1.249 ± 0.043. Thus, the over- 
all model can predict affinity of sets of 5 training set 
ligands using only the information from 12 of the 17 
compounds in the 5HT 7 training set. 

The contour map derived from the CoMFA model 
illustrated several potentially important aspects of 
drug interactions at 5HT 7 receptors. As with, other 
recent CoMFA models for Di and D 2 dopamine recep- 
tors [2-4] steric features near the cationic nitrogen are 
of some importance in differentiating ligand-binding 
affinity for the 5HT 7 receptor. A second impor- 
tant characteristic is that regions in which positive 
charge is favored include the vicinity of the hydrogen- 
bonding nitrogen, and the regions near the live- and 



six-membered rings. Steric features, especially re- 
gions in which steric bulk is favored For high affinity, 
appear prominently in the 5HT 7 model. Consistent 
with what we recently observed for the D2 dopamine 
receptor [4], addition of a hydrophobic component 
slightly worsened the 5HT 7 model. This suggested a 
potentially important role for this factor in positioning 
the alkyl side chain of the highest affinity compounds 
properly in the receptor active site. 

There are a number of papers modeling the inter- 
actions of drugs with various 5HT receptors (see [22| 
for review). Wishart el al. [22] noted that the 1994 
model of Kuipers 137] was able to account for the 
high-affinity of 8-OH-DPAT to the 5HT, A receptor 
as well as its relative selectivity for that receptor vs. 
those of the 5HT 2 subfamily. As expected, the Kuipers 
[37] model favored a three point binding interaction 
between the cationic nitrogen of the drug and the Asp 
of transmembrane helix III. The investigators also sug- 
gested that two hydrogen bonds were formed between 
the drug and a Ser and Thr of helix V. The authors 
further suggested that the ability of one of the bulky 
n-propyl groups of 8-OH-DPAT to form hydrophobic 



XP009032842 



908 

interactions with residues in helices HI and IV of the 
5HT| a receptor that were not accessible to the drug at 
the 5HT2 receptors may have provided part of the basis 
for selectivity. This result implies that the 5HT7 recep- 
tor should also have homologous residues located in 
helices III and IV. Indeed, a preliminary examination 
of the 2D structure of the 5HT7 receptor suggests that 
appropriate residues are well positioned in helices III 
and IV. It will be important to develop a protein ho- 
mology model for the 5HT7 receptor to determine the 
extent to which this is actually possible. 

Based also on the high affinity of 8-OH-DPAT for 
5HT)a and 5HT7 receptors, it may be reasonable to 
extrapolate from some mutagenesis data obtained for 
the SHTia receptor to the 5HT7. It was proposed 
that Asn396 and Ser393 of the rat 5MTi A receptor 
played an important role in the selective affinity of 8- 
OH-DPAT because mutations at either residue reduced 
drug affinity [38]. Accordingly, it is predicted that ho- 
mologous residues in the guinea pig 5HT7 receptor 
would also be exposed to the surface of the binding 
pocket and in a position to influence agonist binding. 
The guinea pig 5HT7 receptor (present study) has an 
Asn in position 389 that seems suitable (homologous 
to the rat Asn 396). The receptor also has three Ser 
residues at positions 360, 378, and 398, which also 
might fulfill the same role as the Ser 393 in the rat 
5HT| A receptor [39, 40]. 

The CoMFA mode! derived here for affinity at the 
5MT7 receptor provided a useful basis for developing 
n simple search query for probing 3D chemical data- 
bases. Limited numbers of hits from each database 
were obtained which were easily rank ordered using 
the similarity index feature present in SYBYL molecu- 
lar spreadsheets. Two compounds with high similarity 
to the template compound and, thus, with significant 
potential as 5HT7 receptor ligands were discovered 
with this procedure. This result offers exciting possi- 
bilities for data-mining new ligands for this important 
receptor. 

It was of interest that the 3D-search/similarity 
index procedures yielded two compounds with 
marked structural similarity to the template R-lisuride. 
Lisuride has high affinity for receptors of the D2 
dopamine subfamily and is a Parkinson's disease drug 
[3, 34, 41]. This suggests that the identified agents- 
may also have some affinity for dopamine receptors. 
The fact that the two potential lead compounds were 
identified as anti-hypertensive agents is intriguing. 
The 5HT7 receptors are expressed in the coronary ar- 
teries and systemic capillary beds as well as in brain 



areas regulating vascular functions fl, 16, 42, 43]. 
Generally, compounds haying the tetracyclic ergoline 
skeleton and derivatives of such compounds possess 
many useful features as potential lead agents, includ- 
ing high affinity and selectivity for the target receptors 
[44]. Compounds such as the two leads identified 
here may have multiple actions on the cardiovascular 
system. 

In conclusion, we have demonstrated the util- 
ity of applying CoMFA to a novel G-protein cou- 
pled serotonin receptor, the 5HT 7 receptor. Similar- 
ities between this CoMFA model and that for the 
Dj dopamine receptor suggest that the CoMFA ap- 
proach can be applied to other systems including D2 
dopamine receptors. It will be of interest to compare 
CoMFA models for 5HT 7( vs. 5HTi a to see if the same 
relationship holds as that observed between the D| and 
Do dopamine receptors. 
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